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Abstract 

The structural and electronic properties of graphene leads its charge carriers to be- 
have like relativistic particles, which is described by a Dirac-like Hamiltonian. Since 
graphene is a monolayer of carbon atoms, the strain due to elastic deformations will 
give rise to so-called 'pseudomagnetic fields (PMF)' in graphene sheet, and that has 
been realized experimentally in strained graphene sample. Here we propose a realistic 
strained graphene bilayer (SGB) device to detect the pseudo-magnetoexcitons (PME) 
in the absence of external magnetic field. The carriers in each graphene layer suffer 
different strong PMFs due to strain engineering, which give rise to Landau quanti- 
zation. The pseudo-Landau levels (PLLs) of electron-hole pair under inhomogeneous 
PMFs in SGB are analytically obtained in the absence of Coulomb interactions. Based 
on the general analytical optical absorption selection rule for PME, we show that the 
optical absorption spectrums can interpret the corresponding formation of Dirac-type 
PME. We also predict that in the presence of inhomogeneous PMFs, the superfluidity- 
normal phase transition temperature of PME is greater than that under homogeneous 
PMFs. 
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Graphene has been an optimal playground to realize the exciting ideas in condensed 
matter physics [l-3] since the carriers in graphene behave like massless Dirac particles. 
Many efforts have been devoted to detecting the striking electronic properties of graphene, 
such as Klein paradox j4j, anomalous quantum Hall effect [5j, and so on. However, it 
was soon realized that the structural and mechanical properties should also be important 
both from theoretical interest and from applications of graphene such as the strain 
engineering in graphene Q-Q. It has been shown recently that specific forms of strain 
oroduce a strong PMF in graphene, which effectively break the time-reversal symmetry 
9]. The strain-induced PMF is expected to produce PLLs, and consequently, the quantum 
Hall effect, even in the absence of external magnetic field 10j. These intriguing properties 
extend to graphene multilayers. Recently, the shear mode in graphene multilayers has been 
observed in Ramman spectroscopy experiment [ijj] , and the high temperature Bose-Einstein 
condensation and superfluidity of indirect excitons or electron-hole pairs in Graphene n- 
p bilayers 16Ml9| are also predicted. Differing from other bilayer n-p systems, such as 
coupled semiconductor quantum wells |20N22l] . n-p type SGB separated by a dielectric layer 
should be an ideal setup to create PME due to the PMF induced by strain rather than the 
applied magnetic filed. Because of the difference of ripples or elastic deformations in each 
graphene layer sample on a substrate, in fact, it is more realistic to fabricate different or 
say inhomogeneous PMFs in electron and hole layers, which is hard to realize with external 
magnetic fields. 

On the basis of these interesting and reasonable ideas, in the present work, we will theo- 
retically analyze the PME properties in SGB device without applying magnetic field. In the 
following, we show that the formation of PME in SGB could be determined from the optical 
absorption selection rule of PME, which is related to the imbalance parameter ^=\/ B^jB e s 
of the strain-induced PMFs B h s and B e s suffered by Dirac holes and electrons in SGB. We also 
find that comparing with homogeneous case (7=1), the Kosterlitz-Thouless (KT) transition 
critical temperature will be improved by a factor of [(2+7+7 -1 ) /4] 2 . Moreover, we suggest 
some technical skills for designing the SGB setup that can make use of this realization to 
detect PME. 

As illustrated in Fig. 1, the system we considered here is that of two parallel strained 
graphene layers separated by an insulating slab of Si02- By varying the chemical potential 
via tuning the bias voltages of the two gates located near the corresponding graphene sheets, 
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FIG. 1: (Color online) Basic Scheme of the proposed SGB device. Two strained graphene 
monolayers are separated by dielectric spacer. Electron/hole carriers in each layer are created by 
n(p)-type doping or applied external gates. Indirect PME can be generated by PMFs in SBG. 
Uniform PMF perpendicular to the graphene surface is created in a circular arc strained graphene 
ribbon. 

the spatially separated electrons and holes in different layers can be generated. In order to 
obtain the analytical PLL expression, we suppose both graphene sheets are bent to circular 
arc [ll|, which can lead to PMFs B e J h = 8/3c e ^ h /a (in unit of H/e=l), and the corresponding 
fictitious gauge field vector potentials can be chosen as A e J h =B e J h (0,x e / h ). Here, do is the 
lattice constant, 0—— <91nt/<91nao~2 and t denotes the nearest-neighbor hopping parameter, 
and c e l h is the numerical constant representing the strength of strain in electron/hole layers. 
In general c e ^c h , and thus B e s ^B^ in SGB device, which results in the fact that 7 departs 
from 1. 

This SGB system is described by the Hamiltonian H=Hq+U (r 6 — r h ), where 

H =hv f [nla 1 ®aQ-nla 2 ®<Jo+irx a Q® (T i-' K yVQ®V2\ , (1) 

with Fermi velocity m-s -1 , ir e / h =p e / h ^f^A e / h } and Pauli matrices <jj (i=0, 1,2). 

U (r e — r h )=— e 2 /t\J\r e — r h \ 2 +d 2 is Coulomb interaction (CI) between the pair of spatially 
separated electron and hole, where d is the spacer thickness and e (~ 4.5 for Si02) denotes the 
dielectric constant of spacer. In the following text, the relative CI strength notation is taken 
as A= (e 2 /e) /(Hvf). Taking some coordinate transformations, the free part of Hamiltonian 
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FIG. 2: (Color online) Dispersions of the first few PLLs for fully occupied case (upper 
panels) and partially occupied (lower panels) case, as functions of P with different 7. 

The dashed horizontal lines denote the non-interacting PLLs. The parameters are chosen as ^=0, 
A=0.49, d=0.2, and R=0. 



H can be written as 



~V2(lJ h f 2 



( l h c + \J_ \ 

l hC \ l e d 

l e c- l h c + 

V Z e c_ l h c\ / 



(2) 



with /, 



e /h=\J h/ eB e J h the pseudomagnetic lengths in electron/hole layer (in following text, 
we set the length unit as y/lJh) and the harmonic lowering operators c± (see Supplementary 
Information part I). 

To obtain eigenvalues E n+ ^ n _ of the total Hamiltonian H, the CI part U needs to be 
treated as a perturbation in the first order o(^+, n'_\U\n + , n_) , where \n + , n_) Q is the eigen- 
states of H corresponding to the eigenvalues, i.e., the PLLs of H without CI 



=— 

n + ,n. ^ 



_ v ^(l+ 7 - 1 )- S -v^M(l+7) 



(3) 



in units of Hvj/y/lJh, with s.|-=sgn(ri-i-). The imbalance parameter 7 will reduce the degen- 



eracy of PLLs of H . For instance, |1, 0) and |0, — 1) are degenerate when 7=1, while these 
two states turn to be nondegenerate when 77^1. The PLLs E n+:n _ of PMEs are not only 
dependent on 7, but also related to the effective momentum P=2^!±i^. an( j coordinate 
R= lhr , e tr of the mass center of PME since U is related to P and R after transforma- 
tions (see Supplementary Information part I). The location of the chemical potential will 
determine the attainable possible PLL indices for electrons and holes. For convenience, in 
this work we choose the notation /i describing the highest filled PLL and consider two cases: 
(i) the electron-PLLs with n + >fi are unoccupied and the hole-PLLs with n^^fi are fully 
occupied; and (ii) the electron-PLLs with n + >fi are unoccupied and the hole-PLLs with 
n_</i are fully occupied, while the PLL at [i is partially occupied. 

The spectrums of the lowest five PLLs as a function of P are demonstrated in Fig. 2. 
Figures 2a— 2e correspond to the fully occupied cases while Figs. 2f— 2j are for the partially 
occupied cases. The PME dispersion strongly depends on the imbalance parameter 7 (see 
Fig. 2). On one hand, if the difference in the elastic deformations between electron and 
hole layers is large, i.e., 7 departs from 7=1 greatly, the PME dispersion will accelerate 
the transition process from density wave forms to particle-hole pair forms with increasing 
the momentum P. On the other hand, when 7 takes a moderate value around 1, the 
PME dispersion behaves like that of magnetoexcitons in perfect graphene bilayers under an 
external magnetic field. Besides, we also found that, taking the partially occupied cases 
as an example, at much larger 7 the lowest PLLs behave like electron type |n + ,0) (for 
example, Fig. 2j) while at much smaller 7 the lowest PLLs behave like hole type |0, n_) (for 
example, Fig. 2f). Differing from the magnetoexcitons in an external magnetic field, the 
PME dispersion is also related on R. By increasing R, the PME dispersion at smaller P 
turns to become a particle-hole pair form rather than a density-wave form (Supplementary 
Figure SI). 

The optical absorption spectroscopy analysis is an instructive method in studying the 
properties of magnetoexcitons, because it can determine the knowledge of magnetoexci- 
ton's energies and wave functions at P=0. The particle-hole excitation energy should be 
revealed directly by the resonant peaks in the magneto-optical absorption spectroscopy 
R a b(R, w)oc J2 n+ n- I ( n +> n - I^xIa*) A 4 ) f ■ I* 1 the absence of inter-PLL CI which is much weaker 
than intra-PLL CI, for linear polarized light A=A~k, the optical absorption selection rule 
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FIG. 3: (Color online) The optical absorption curves of PMEs with //=0. a, Fully and b 

partially occupied cases as functions of photon energy u (h=l). The spacer thickness is set to 
d=0.2. A=0.49. 

for PMEs is analytically expressed as 

(n+, n-\v x \fi, //)=ci5| n+ [ ) [ /Lt |5[ T j_| i | jU |_i+c 2 5| n+ | i | jU |_i5| n _[ ) | /Lt | 

where v x =dHo/frdp x is the velocity operator for free part of Hamiltonian, and the 
coefficients ci= -(l+7)s_(s (U s + +l)C, c 2 =(l+7 _1 )(s M +s_)C, c 3 = — (l+7)(s /x +s+)C, and 
c 4 =(l+7- 1 )s + (s M s_ + l)C with C=v f (V2) Sn+ ' 0+Sn -' 0+2SlJ "°~ 6 (see Supplementary Information 
part II). This optical absorption selection rule formula fll]) offers a better way to analyze 
PME absorption resonant peaks in spectroscopy of SGB device. Notice that the optical 
absorption formula fll]) also depends on R since the PME spectrum is related to R. (For 
briefness, we just consider the R=0 case in following text.) 

Especially, for the fully occupied case with /x^O, the selection rule can be further simplified 

as 

(n + ,n_\v x \fi,fi) = c 4 5 n+tll+1 5 n _^, (5) 

which indicates that there is only one resonance peak occurring in optical absorption spec- 
trum, which corresponds to the formation of PME. For example, the selection rule for fully 
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occupied case of /i=0 with the ground state 1 0, 0) promises that just one resonance peak 
that corresponds to the formation of PME state |1, 0) occurs, see Fig. 3a. From Fig. 3a one 
can also clearly observe that with the imbalance parameter 7 decreasing (increasing) from 
1, the resonance peak moves towards high (low) frequency region, see the red (blue) curve 
in Fig. 3a. 

This effect of inhomogeneity of PMFs can also be seen for the partially occupied cases, 
which is illustrated in Fig. 3b. According to the general selection rule fll]), the nonzero 
elements of (n + , ri-\v x \0, 0) only occurs at two PMEs states |1,0) and |0, — 1), and thus in 
this case, there are two resonance peaks occurring in Fig. 3b. Interestingly, with decreasing 
(increasing) 7 from 1, the resonance peak of |0, — 1) moves towards low (high) frequency 
region, in contrast, the resonance peak of |1,0) moves towards high (low) frequency region. 
The results of cases with /x^O are similar to these of /i=0, and an example for is shown 
in Supplementary Figure S2. 

At this stage, we have not considered the effect of the inter-PLL CIs on the optical 
absorption of the PMEs. Because the inter-PLL CI mixes the noninteracting PLLs, the 
absorption phenomenon should appear at the energy of every PME states. Consequently, 
additional resonance peaks, except of the main peaks without the inter-PLL CI, will appear 
in the optical absorption spectrum (not shown here). The magnitudes of these additional 
resonant peaks are very tiny and tend to disappear by increasing the spacer thickness d. 
Therefore, based on the optical absorption selection rule analyzed above, one can see that 
the optical techniques should allow detection of PME locally as well as the CI effects in the 
present SGB device. 

PMEs are also ideal objects in exploring the Bose-Einstein condensation since they behave 
as neutral bosons at low densities. Motivated by its importance both from basic point of 
interest and from application of graphene-based electronics, now we turn to address this 
issue by presenting an attempt at the theoretical evaluation of effective mass the effective 
mass mg s (B s =yJ B^B^ herein) and superfluid-normal state, i.e., KT transition temperature 
T c of PME in this SGB system. As an example, let us now consider the case of PME on the 
(1, 1)-PLL, with the PME energy at small magnetic momenta E^iiP) —^b s (^) + 2mB (cQ ( see 
Supplementary Information part III and ref. 26). In the special limit of d^y/IJh, one can 
explicitly obtain — — fs> which is same as that of magnetoexcitons in external magnetic 
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FIG. 4: (Color online) Calculated KT temperature versus the PMF. The density of PME 
is set as n=4.0 x 10 11 cm -2 and the spacer width d=30 nm. 



field since (d) is independent of 7. Whereas the effective mass of PME should be 

m B = [4/(2+7+7-)] 2 



(6) 



where the prefactor [4/ (2+7+7 -1 )] 2 will reduce to unity in the case of homogeneous PMFs 
(i.e., 7=1), thus that the effective mass of PME in SBG is [4/ (2+7+7- 1 )] 2 times of that 
of magnetoexcitons in perfect bilayer graphene under an external magnetic field. It is this 
difference between the effective masses that results in the difference in the KT transition 
temperature T c , which reads 

1/3 

(7) 




32 



-[(sm Bs k B T«) I {TtVn)f +1+1 



QO 

-[(sm Bs k B T0)/(^n)] 3 +1-1 



1/3' 



2^0 
2i73' 



where the auxiliary quantity is T 6 °=-^ ( 6s((3)m B 



1/3 



with fi - 



wh 2 r 



srriB s log sS 4 e 2 /(27rnm 2 « 4 



being 



the chemical potential of system, s=4 the spin degeneracy of PMEs in SGB, and ((3) the 
Riemann zeta function. Here, n = n s + n n with n n the normal component density, and 
n s the superfluid density. Because fio approximately inverse of the effective mass, a rough 
consequence can be immediately obtained that the critical temperature of KT transition of 
PME in SGB will be greater [(2+7+7 -1 ) /4] 2 times than that of magnetoexciton in bilayer 



graphene in an external magnetic field 18] . Typical results are exhibited in Fig. 4. Thereby, 
one can conclude that increasing the imbalance of the strains (7 departs from 1 sufficiently) 
in different layers in SGB device may be one effective method to promote the KT transition 
temperature of PME. 

In order to design the required SGB device in experiment, firstly, one should deposit 
graphene ribbons onto a elastic substrate and then deform it by bending it into a circular 
arc automatically (Fig. 1). Following that, we can transfer it onto the Si02 thin film and 
then throw off the elastic substrate by heating or dissolving it via chemical agent carefully. 
A uniform PMF could be created in graphene ribbons deposited on SiC>2 thin film since the 
strain distribution in the substrate would project onto graphene ribbon. In reality, the PMFs 
arisen from strain in upper and lower bent graphene ribbons, however, may be different be- 
cause of different automatic deformations of them. This is why we discuss the importance 
of the imbalance parameter 7 here. After that one can grow or evaporate metal electrode on 
the other side of SiC>2 thin film, and finally package the device for measurements. The ex- 
perimental capability to produce high quality graphene sample and the feasibility to create 
strong PMF in graphene, together with the technical advances for designing graphene multi- 
layer devices make it possible to perform very original optical and superfluidity-normal states 
KT transition experiments of PME in graphene without external magnetic field. Moreover, 
the SGB device proposed here could also be used to detect the (valley polarized) uncon- 
ventional fractional quantum Hall effect of charged Dirac-type electron-hole fluid or a Bose 
condensate of PME in the absence of external magnetic field by measuring transport coeffi- 
cients and that, since fixed PME number corresponds to fixed valley polarization in the SGB 
device, transport studies of the SGB system will provide unique information which cannot 
be observed by direct study of conventional two-component electron-electron/electron- hole 



systems 



23]. 



In summary, in this work we propose a realistic SGB device, and based on the physics 
analyzed here, this setup can be made to detect the PME which is associated with the elastic 
deformations or curvature in graphene. The cases discussed in this paper are only illustrative 
examples of the feasible setup and applications of the SGB device, and the technical skills 
mentioned above are achievable in current experimental capabilities. Therefore, we hope 
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these predictions and the suggested SGB device can be realized sooner rather than later. 
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I. Let us explain first how to derive the Hamiltonian (2) of the main text. Taking the 
coordinate transformations 

p= JeP e +lhP h p= l eP e -l hP h 
h+lh h+lh 

„ l h r e +l e r h J h r e -l e r h 

R= ' ,r=2 h e , (S2) 



we can rewrite h e J h in Eq. (1) as 



h o /h= { K e/h ( Y ±Px ) +i 



(S3) 



where K e / h — (l e +h) I fae/h), and £=1/ (2lJ h ). Here, r y—l e /l h —-^/B~f]B~l describes the inho- 
mogeneity of PMFs suffered by electrons and holes in different layers. In the limit of 7=1, the 



problem of PMEs in SGB is reduced to that of PMEs in a uniform PM 



with those on the magnetoexcitons in an external magnetic field 



which are consistent 
2]. And then one can 



introduce the transformation S=e lXy ' lelh , for which P x S=P x +y/lJh, S j; PyS=Py+X/l e l h , 
and shift r^-r—l e lhZxP, h e J h can be further written as 

hf=K elh {{± Px +Cy) ±i (PyTtx)} . (S4) 

The relative coordinate of electron and hole in the x-y plane r e —r h becomes 
f=\ (7— 7 _1 ) R+7 (2+7+7 -1 ) (r—lJhZxP). Finally we define the harmonic lowering oper- 
ators 

a=VUhPx~i fi^ ,b=\flJhPy-i /rr , (S5) 



and their combination 

c± = ±(a±ib)/V2. (S6) 



12 



Substituting them into Eq. (1), we can obtain the Hamiltonian (2) of the main text. The 
eigenvalues, i.e., the PLLs of H without CI are given by Eq. (3) and the corresponding 
eigenvectors are given by 

( s + s_$| n+ |_i,| n _|_i(r) \ 
s_$| n+ |,| n _|_i(r) 

V ^|n + |,|n_|(r) 



\n+,n.) =2 r ' 



(S7) 



/ 



where Tf= Sn+ ' 0+S " '° 2 , s ± =sgn(n ± ), and $ nil „ 2 (r) 



\'z\ 



n jriV rM Ll^)e- r l with 



lz—ni—ri2, n-— min(ni, 712), S (l z )—sgn(l z ) lz — >1 for l z —0, and L(x) the Laguerre polyno- 
mial. To obtain eigenvalues of the total Hamiltonian H, we need to solve the following 
equation: 



0=det 



5 n+:n , + 5 n _ :n ,_(EM n _-E) 



+o(n' + ,n'_\U (?) |ra+,ra_) | 



(S8) 



Here, the intra-PLL component of the CI is defined as o(n + ,n^\U |n + ,n_) , while the inter- 
PLL component is defined as o{n' + , n'_\U \n + , ra_) , where \n' + ,n'_} 7^ \n + ,n_) . The PLLs 
as a function of R are presented in supplementary Figure SI. 

II. From the Fermi's golden rule, the optical absorption for photons of frequency u> yields 



2tt 



R ab (R,ou)=—J2\( a \- A - v \U) 5(e ( «\R)-e -hLu), 



(S9) 



where |a) is the set of quantum numbers describing a particle-hole excitation and \0)=\/i, /i) 
is the ground state in absence of particle-hole excitations with £ being the corresponding 
ground state energy. v=dH /Hdp is the velocity operator for free part of Hamiltonian, and 
A is the vector potential. For linear polarized light A—Ak, choosing a Lorentzian type 
broadening, one has 

r/2 



Rab(R, w)aj] I (a|u x |0) | 



eS?(-R)-£o-M 2 + r2 / 4 ' 



(SIO) 



where T denotes the linewidth, and | (o;|f ^ j 0) | 2 = J2n + >^ En-^ / ^ r * s f° r ^ un occupation 
while Yln + ^n E n _^ / Tdr is for partial occupation. The parameter T is defined as 



(Sll) 
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Supplementary Figure SI: Dispersions of the first few PLLs as functions of P with 
different R. The imbalance parameter 7=0.5 for a-e while 7=1.5 for f-j, and other parameters 
are same as Figure 2 in main text. 
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Supplementary Figure S2: The optical absorption curves of PMEs with fj,—l. a, Fully 
and b partially occupied cases as functions of photon energy uj (h=l). The spacer thickness is set 
to d=0.2. /3=0.49. The inset in b is the enlarged dispersion curves between w=0.8 and w=1.2. 



where C" + n _ is the projection of PME state \a) on the basis state \n + ,n_). In the absence 
of inter-PLL CI which is much weaker than intra-PLL CI, the PME state \a) is one special 
basis state and C% n =^a>,(n + ,n_)> same as that in the noninteracting case. As a result, we 
can analytically obtain the optical absorption selection rule for PMEs, Eq. (4) in main text. 

As an example, we additionally plot the optical absorption spectrum with fully and 
partially occupied cases in supplementary Figure S2. The selection rule for fully occupied 
case with /i=l promises that there is only one resonant peak to occur, which corresponds to 
the formation of PME state |2, 1) by absorbing a photon quanta ?vjj. For partially occupied 
case, however, there are three non-zero elements according to selection rule (4), which are 
|1,0), 1 1,-2), and |2, 1). Therefore, there are three resonant peaks appearing in optical 
absorption spectrum, see supplementary Figure S2b. Similar to those in Fig. 3 in main 
text for the fi=0 cases, the imbalance of strain-induced PMFs changes the location of the 
resonant peak. When the imbalance parameter 77^1, these resonant peak moves in a more 
complex way. 
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III. In the limit of d-^-oo and B s — >oo, the PME energy can be approximated by only 
considering its zeroth order energy part En,m- However, if the PMF is about 10~20T, CI 
can be treated as a perturbation since CI e 2 /ed is just several times less than hvf/y/lJh- In 
the absence of inter-PLL CI, by substituting the approximate relation 

(( nmP Mf)\nmP))=eSl+ m ^ (S12) 

for R=0 into Eq. (S8), we can get the asymptotic dispersion law of a PME for small 
magnetic momenta. Here the notation ((nmP\U(r)\nmP)) denotes the averaging by the 
two-dimensional harmonic oscillator eigenf unctions $„ jm (r), and Enm is the binding energy, 
and finally one can obtain the Eqs. (6) and (7) in main text. 
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